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HOMOGENIZATION OF GENERALIZED SECOND-ORDER ELLIPTIC DIFFERENCE 

OPERATORS 


ALEXANDRE B. SIMAS AND FABIO J. VALENTIM 

Abstract. Fix a function W{xi, ... ,X(i) = where each Wk : M ^ M is a strictly increasing 

right continuous function with left limits. For a diagonal matrix function A, let VAVvi^ = ) 

be a generalized second-order differential operator. We are interested in studying the homogenization of 
generalized second-order difference operators, that is, we are interested in the convergence of the solution of 
the equation 

\UN - 

to the solution of the equation 

Xu — VAVwu = /, 

where the superscript N stands for some sort of discretization. In the continuous case we study the problem 
in the context of W-Sobolev spaces, whereas in the discrete case the theory is developed here. The main 
result is a homogenization result. Under minor assumptions regarding weak convergence and ellipticity of 
these matrices , we show that every such sequence admits a homogenization. We provide two examples 
of matrix functions verifying these assumptions: The first one consists to fix a matrix function A with 
some minor regularity, and take A^ to be a convenient discretization. The second one consists on the case 
where A^ represents a random environment associated to an ergodic group, which we then show that the 
homogenized matrix A does not depend on the realization uj of the environment. Finally, we apply this result 
in probability theory. More precisely, we prove a hydrodynamic limit result for some gradient processes. 


1. Introduction 

In the ’50s William Feller introduced a more general concept of differential operators, that is, operators 
of the type {d/dW){d/dV) where, typically, W and V are strictly increasing functions with V (but not 
necessarily W) being continuous. In this paper we are interested in the formal adjoint of {d/dW){d/dV), 
which is simply {d/dV){d/dW)^ in the case V{x) = x is the identity function. For more details on Feller’s 
operators, we refer the reader to 011 nn]. 

Recently, the formal adjoint operator (d/dx){d/dW) and some non-linear versions were obtained as scaling 
limits of interacting particle systems in inhomogeneous media. They may model diffusions with permeable 
membranes at the points of the discontinuities of W, see [I10[710[I4] for further details. In [14], for instance, 
the author introduces an extension of the formal adjoint operator to higher dimensions and provides some 
results regarding this extension. Informally, fix a strictly increasing right continuous functions with left 
limits and periodic increments, Wk : R —>■ R, fc = 1,..., d, and let W{x) = for x G R'^. 

The generalized Laplacian operator is given by Lw = Note that Wk{x) = x, k = 1,... ,d, 

corresponds to the classic Laplacian operator. Furthermore, In m the authors studied the space of functions 
/ having weak generalized gradients 'Vwf = which they called the W-Sobolev spaee. 

Several properties, that are analogous to those in classical Sobolev spaces, are also obtained. Recently, in 
m, W-Sobolev spaces of higher order were defined and results on elliptic regularity were obtained. 

Our goal in this paper is to approximate a solution of a certain continuous partial differential equation 
driven by a generalization of the laplacian by considering a sequence of discrete functions given by solutions 
of discrete versions of such equations. 

A first contribution of the current paper, presented in Section 14.11 is to build a theory in discrete setup 
that is analogous to the theory introduced in [12]. In particular, we introduce the VF-interpolation and 
present some results in order to get connections between the discrete and continuous Sobolev spaces. 
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The main contribution are the results related to homogenization presented in Section [S] Informally 
speaking, we considered discrete problems of the form Am at — where are sequence of 

diagonal matrices functions, then we define a matrix A as being a homogenization of the sequence of matrices 
A^ if some convergences, in appropriate spaces hold, and also if the limit is related a continuous problem 
associated the matrix A. The main result of this paper, namely Theorem IS.Sl can be summarized in following: 
under minor assumptions regarding weak convergence and ellipticity on these matrices A^, we show that 
every such sequence admits homogenization. 

This result can be incorporated in several situations. We provide two examples that cover many of them. 
The first consists in fixing a matrix function A with some minor regularity, and let A^ be a convenient 
discretization. Theorem 15.81 ensures the convergence of the solutions un to uq, which is associated the 
problem with A. The other example consists of an example of stochastic homogenization. Indeed, let A^ be 
a sequence of random matrix functions satisfying an ergodicity condition. Then, we show that, this sequence 
admits homogenization, and that the homogenized matrix A does not depend on the realization uj of the 
random environment. The focus of this approach is to study the asymptotic behavior of the coefficients for 
a family of random difference schemes, whose coefficients can be obtained by the discretization of random 
high-contrast lattice structures. In this sense, we want to extend the theory of homogenization of random 
difference operators developed in m, which was also tackled in seminal paper by Kozlov [5]. We also want 
to generalize its main Theorem (Theorem 2.16) to the context in which we have weak generalized derivatives 
in the sense of IT-Sobolev spaces. 

Finally, as an application of this theory, we prove a hydrodynamic limit result, which is, per se, an 
interesting result, since it extends the main result obtained by M- More precisely, one must observe that 
this is a non-trivial extension, since the fact that the matrix A can be non-constant creates a large difficulty 
to the method presented there. More precisely, the method consisted in studying the spectral decomposition, 
and the resolvent operator of the operator Cw = X]fc=i ^x^dwh- Thus, in our setup, that allows the matrix 
function A to be non-constant, would force one to look for spectral decomposition and the resolvent operator 
of the operator J2k=i dxk{o,kdwk)- 

The remaining of the article is organized as follows: in Section [2] we provide a brief review on W-Sobolev 
spaces; in Section [3l we define a new space of test functions needed in subsequent sections; in Section |4] we 
introduce a discrete analogous to the continuous W-Sobolev spaces, and then, interpolation and projection 
results connecting the discrete and continuous versions; in Section [5] we provide main results of this article; 
in Section we apply the results of this article to prove a hydrodynamic limit for gradient processes with 
conductances in random environments. Finally, we provide an Appendix with a proof to an auxiliary result. 


2. W-SOBOLEV SPACE 

In this Section we recall some notation and results of [3 m m]. Denote by = (R/Z)'^ = [0, l)'^ the 
d-dimensional torus and fix a function W : —)■ R such that 


( 1 ) 


d 

W{xi, ...,Xd)=^ Wk{xk), 

k^l 


where each Wk : M —^ R is a strictly increasing right continuous function with left limits (cadlag), with 
periodic increments, in the sense that for all u S R, Wk{u + 1) — Wk{u) = VFfc(l) — Wfe(O). 

Define the generalized derivative dw^ of a function / : —>■ R by 


( 2 ) 


dwkfixi,.. .,Xk, ...,Xd) = lim 

e^O 


fjxi,.. .,Xk + c, ■ ■ ■ ,Xd) - fjxi,. . .,Xk,. ■ ■,Xd) 
Wkixk + e) - ITfc(xfc) 


if the above limit exists. Denote the generalized gradient of / by Vwf = {dwi 9wdf ); if the generalized 
derivatives dw^ exist for all fc = 1,..., d. 

Let us remember the definition of the space idi,vv(T‘^), called VF-Sobolev space. We denote by (•, •) the 
inner product of the Hilbert space and by || • || its norm. Let be the Hilbert space of 
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measurable functions / : T'^ —>■ K such that 




Jfd 

where Wk) represents the product measure in obtained from Lesbegue’s measure in ^ and the 

measure induced by Wk in T: 

d{x^® Wk) = dxi ■ ■ ■ dxk-i dWk dxk+i ■ ■ ■ dxd- 


w '' ft,y X ^“^ft,- 

Denote by (•, ■)x’‘^Wk the inner product of 


{f,9)x>^®Wk = I fix)g{x) d{x^(g>Wk), 

JTd 

and by || • lU^igiu/;, the norm induced by this inner product. 

The set Awk of the eigenvectors of forms a complete orthonormal system in L^(T), (see [5]). Let 

d 

•^W = {f ^ f{xi, ■ ■ ■ ,Xd) = Y\_ fk{xk), fk & Awk}j 

fc=l 


and denote by Dw := span{Aw)- Define the operator £w ■ Dw := span{Aw) as follows: for 

/ = nfc=i fk G Aw, 

d d 

(3) Cw{f){xi,.. .Xd) = E n fj {Xj)Cwk fkixk), 

k=l j = l j'/fe 

and extend to Dw by linearity. In particular, Dw is dense in L^(T'^) and, the set Aw forms a complete, 
orthonormal, countable system of eigenvectors for the operator Cw- More details can be found in [Ill- 
Denote by closed subspace of L'^k^-^^^{T‘^) consisting of the functions / that have 

zero mean with respect to the measure d{x^® Wk)'- 

[ f d{x^® Wk) = 0. 

Jfd 

The function g € has W-generalized weak derivative if for each k = 1,..., d there exists a function 

Gk € L'^k^Wk satisfying the following integral by parts identity 

(4) f (dxkdwK f) 9 dx = - [ (dwkf) Gkd(x'^<g>Wk), 

Jjd Jjd 

for every function / G Dry. 

Denote by iJi^vy the ID-Sobolev space, the set of functions in having IT-generalized weak 

derivative. Note that Dw C Hi^wO^‘^)- Moreover, if g G Dw then Gk = dwk9- For this reason for each 
function g G Hi w "we denote Gk simply by dwk 9 , we call it the fc-th generalized weak derivative of the 
function g with respect to IT. 

In [12] it is shown that Gk is unique and the set iLi,vy(T'^) is a Hilbert space with respect to the inner 
product 

d ^ 

if, 9 ) 1 ,w = (/,5)+E / i 9 wkf){dwk 9 ) d{x'"(S)Wk). 

fe=ldT<i 

Furttisrinorc, let (^TT ) be the duRl sp&cc to (T'^), that is, is the set of bounded linear 

functionals on Hi_wO^'^)- The following characterization of H^(T‘^) can be found in [12] : 

Lemma 2.1. / G H^{T‘^) if and only if there exist functions fo G L^(T'^), and fk G L'^k^Wk 
that 


d 

k^l 
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( 5 ) 


in the sense that for v € Hi ,wir^) 


Furthermore, 




/(^) := [ fovdx + Y] [ fk{dwkV)d{x^®Wk). 
inf I l/ol^dx + |/fc|^d(a;''(g) VFfc)^ ; / sate/ies dS]) | . 


3. A NEW SPACE OF TEST FUNCTIONS 

We will now define a new space of test functions that will be needed in the proof of the hydrodynamic 
limit in Section [51 and also to prove the Compensated Compactness Theorem (Lemma 15.21) . 

In this paper we deal with discrete approximations of functions defined on the torus T^. To obtain 
a hydrodynamic limit, we will eventually need to apply these discretizations to test functions of the W- 
Sobolev space. The problem is that these test functions are defined as functions in and thus, the 

discretization procedure obtained by restricting the test function to is not well-defined. Therefore, we 
will define a new space of test functions suitable for discretizations, which we will denote by Note 

that in [12| we considered another space for test functions, namely Dw, whose definition was given in the 
previous Section. 

Another motivation for such definition is of independent interest: to develop a classical theory of functions 
that admit generalized derivatives. In fact, these test functions will be such that one may apply the operator 
dxkdwk in the classical sense. 

To this end, recall, from o, the definition of W. For / : T —>• R, let D{f) be the set of its discontinuity 
points. For k = I, ..., d let also Cw^ (T) be the set of cadlag functions / : T -> R such that D{f) C D[Wk)- 
We endow Cw^^) with the sup norm || • ||oo. Equation dJ) in the one-dimensional case becomes 

rf/ / s, f{x + e)-f{x) 

-(x) = lim ■ 


dWk^""' <i^oWk[x + e)-Wk[x)' 
if the above limit exists. Let be the set of functions / in such that is well-defined and 

differentiable, and that belongs to C'wfc(T). 

In [5], it is proved that Dwk is the set of functions / in Cwki'^) such that 

(6) f{x) = a + bWk{x) + [ dWk{y) f g{z)dz 

J{0,x] Jo 

for some function g in C'wic(T), where a, b are real numbers satisfying 

(7) bWkil) + [ dWk{y) [ g{z)dz = 0, / g{z)dz = 0. 

JT Jo Jt 

The first requirement corresponds to the boundary condition /(I) = /(O) and the second one to the boundary 
condition df /dWk{l) = df /dWk{0). 

Let us now define a d-dimensional counterpart to the sets C'wfc(T) and respectively: 

C'iu(T‘') =span{/;/ = /i(g)---(8)/d, /fc G (T), fc = 1,..., d} , and. 


£>w(T'^) = span{/;/ = /i 


' fd, /fc G Dwfcjfc = 1, ■ • ■ , 


( 8 ) 

( 9 ) 

where ® denotes the tensor product. That is, for functions /i G Dwn ■ ■ ■, fd & we have / = /i®- • -^fd 

is such that 

d 

f{xi ,... ,a:d) = fk{xk), 

k^l 

and span(A) stands for the linear space generated by the set A. In other words, Cn/(T^) and are 

the spaces generated by functions of the form fi ^ ‘ ^ fd- 
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Note that is “almost” the tensor product of the spaces Dwi, ■ ■ ■ ,'^Wd: the reason why it is not, 

is that the tensor product of topological spaces is defined as a closure of the previous space, and we do not 
want to take this closure. 

Remark 3.1. Note that if f G £)w(T'^), then admits a partial derivative in the kth direction, therefore, 

if Qk G C\y(T‘^) is such that its partial derivative in the kth direction exists, then, the function akdyy,,f also 
admits a partial derivative in the kth direction and we have 

dxk (akdwkf) = {9wkf){9xkO,k) + akd^^^dwkf 

in the strong sense. 

The above remark motivates the following definition. Let M 14 / C be the space of functions 

a = (oi,..., Od) such that for every k = 1,..., d, Ok admits a partial derivative in the /cth direction. By 
convenience, we will also say that a diagonal matrix A = (okk) belongs to Mw to mean that the function 
a = (oii,..., Odd) belongs to Mw- 

Remark 3.2. It is easy to see that if f G one may apply the operator dx^dwk ihe classical 

sense. Actually, we have that for f G S)u/(T'^), and a G Mw, VAVw/ G C'w(T'^), where 'S/AVwf = 
(akdwkf)- 

The following result shows that can be used as a space for test functions: 

Proposition 3.3. The space SvvCIT'^) satisfies 

i) Sw(ir‘^) is dense in in the L'^-norm; 

ii) 2 )^( 1 ’’^) is dense in C(T'^), the space of continuous functions in in the sup norm || • ||oo. 

Proof. We begin by noting that the first statement follows directly from the second statement. We will, 
thus, prove the second statement. 

Begin by noticing that from Proposition 17.11 in Appendix, SiVt is dense in C'(T) in the sup norm. To 
extend this result to the d-dimensional case, note that, since T is a compact Hausdorff space, we can use 
Stone-Weierstrass Theorem to conclude that 

d 

k=l 

where {^^^1 C'(T) is the tensor product of the spaces C'(T), that is, C'(T) is the closure of the space 

generated by functions of the form /i (8) • • • 0 /d, with fk G C{T),k = 1,... ,d. 

Now, it is easy to see that, since Dw^ is dense in C'(T), is dense in C'(T'^). 

□ 

Let us define a function / : T'^ —R as right-continuous if for all x S T'^, we have lim^j^j, f{y) = f{x), where 
y f X means that for all k = 1,... ,d, pk i Xk- Observe that for G Cwk we have \\-aiy,^^xk fkiVk) = fk{xk). 
Thus, it is clear that for / = /i 0 • • • 0 /d, we have limyj,^; f{y) = f{x). We then have the following remark: 

Remark 3.4. Every function f G C'w(T'^) is right-continuous. In particular, since C C'w(T'^), 

every function f G is right-continuous. 

Lemma 3.5. Let u G C'w(T'^), then u is bounded. In particular, if u G S)w(T‘^), then u is bounded. 

Proof. From the definition of Cvi/(T‘^), there exist bj G M, and functions G CwkO^)j k = l,...,d, 
j = 1,..., m, for some m G N, such that 

m 

u{x) = ^ bju{ 0 • • • 0 u^di^)- 
t=i 

Therefore, it is enough to show that each function is bounded in T. Suppose that u^. is not bounded, 
then, there exists a sequence (t„) in T such that 

\ui{t^)\-^oo. 
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By compactness of T and using the identification T = [0,1), we can find a monotone subsequence such 
that lim^ G T. If tm is increasing, then the existence of left limits of shows that |n^(tni)| bounded, 
whereas if is decreasing, the right-continuity shows that \ is bounded. This contradiction concludes 

the proof. □ 


4. The discrete 1T-Sobolev space 

In this Section, we introduce the notion of discrete IT-Sobolev space. As in the continuous case considered 
in |12j . we obtain similar results for the discrete setup. Furthermore, we introduce the W- interpolation, and 
we also deal with other types of interpolations. Finally, we present some results in order to get connections 
between the discrete and continuous Sobolev spaces. 


4.1. The space Hi ^wO^n)- Denote by = {X/N'LY = {0, ..., — I}'* the d-dimensional discrete torus 

with N'^ points. We will now define and obtain some results on a discrete version of the W-Sobolev space. 

For / : > K. , consider the following difference operators: 3,^, which is the standard difference 

operator on the fcth canonical direction, and is given by 


( 10 ) 



= N 


/ 


x + ek'\ 

N J 



for X e T^; 


and , which is the WT--difference operator given by 


d 


N 

Wk 



w{ag^)-w{j,Y 


for X G T^. 


Denote by L^(T^), and the Hilbert spaces of the functions defined on obtained 

with respect to the following inner products: 


(f,3}N ■= 

{f.9)w,,N:=j^ fY/N)g{x/N){Wi{x + ek)/N)-Wix/N)), 

xeT^ 

d 

{f^9)l,W,N ■={f^9)N -h '^{dwkf^9wk9)wk,N, 

k^l 


respectively. Consider the following notation for their induced norms: 

II/IIl 2 (t^) = {f,f)N, = {fJ)wk,N and ll/ll^i = {fJ)i,w,N, 

respectively. 

These norms are natural discretizations of the norms considered in continuous case. Using the same 
arguments used in |12| to prove Poincare inequality, one can easily prove its discrete version: 


Lemma 4.1 (Discrete Poincare Inequality). There exists a finite constant C such that 








L^T%) 


for all f : jfTff —>■ K, where 

(T%)- 

k=l 

We will now provide results on discrete elliptic equations. We remark that the proofs of these lemmas are 
identical to the ones in the continuous case, and the reader is then referred to [12]. 
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Let A > 0 and A = (akk(x))dxd, x G T'^, be a diagonal matrix function satisfying the ellipticity condition: 
there exists a constant 0 > 0 such that < akk(x) < 9, for every x gT‘^ and every k = 1,. .. ,d. We are 
interested in studying the problem 

(11) = /, 

where u : —?► R is the unknown function, / : N —>■ R is given, and denotes the discrete 

generalized elliptic operator 

( 12 ) Tj^u := Xu-V^AV^u, 

with 

d 

V^AS/^u := Y.d^^lakkixmd^^u). 

k=l 

The bilinear form B^[-, ■] associated with the elliptic operator Tj^ is given by 

B^[u,v] = X{u,v)n + 

(13) 

+ 1^ T,k=i akk{x/N){d^^u){d^^v)[Wkiixk + l)/N) - Wkixk/N)], 

where m, v : t R- 

A function u : —>■ R is said to be a weak solution of the equation T^u = / if 

B^[u,v] = {fTv)]s[ for all r; : R. 

Denote by the subspace of formed by functions / : —>■ R which are orthogonal 

to the constant functions with respect to the inner product (•, ■)i^w,n- Note that i^J‘-^y(T^) is the set of 
functions / : —>■ R such that 

^ E /(^/^) = o- 

Remark 4.2. As seen in Remark \4-.5\ below, a function u is a weak solution of the discrete problem dill) if, 
and only if, it is a strong solution. We chose to present weak solutions of discrete problems because it is easy 
to obtain existence and uniqueness results following this approach, and the results are stated in a very similar 
fashion to those in the continuous case, thus making the analogy with the continuous case more transparent. 

Lemma 4.3. Given a function f : —>• R, the equation 

V^AV^u = /, 

has weak solution u : —>• R */ and only if f G iLj^^(T^). In this case we have uniqueness of the 

solution disregarding addition by constants. Moreover, if u G (K) we have the bound 

where C > 0 does not depend on f nor on N. 

Lemma 4.4. Let A > 0 and f : —)■ R. There exists a unique weak solution u : —>■ R o/ the 

equation 

(14) \u-y^Ay^u = f. 

Moreover, 

where C > 0 does not depend neither on f nor on N. 

Remark 4.5. Note that in the discrete Sobolev space Hi w(T^) we have a “Dirac measure” concentrated 
in a point x as a function: the function that takes value N‘^ in x and zero elsewhere. Therefore, we may 
integrate these weak solutions with respect to this function to obtain that every weak solution is, in fact, a 
strong solution. 
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Remark 4.6. We had to consider diagonal matrices, because a non-diagonal matrix implies the usage of 
an operator of the form dx^dw^, with i ^ j. To apply such an operator in the strong sense, the function Wj 
cannot be discontinuous at any point, and thus this operator loses its main physics motivation, which is to 
create permeable membranes at discontinuity points. 


Motivated by Lemma 12.11 we are able to obtain a characterization of the elements of the discrete dual 
Sobolev space. 

Lemma 4.7. Let f £ the discrete dual Sobolev space. Then, there exist unique functions 

fk ■ —7- R, k = 0,1,... ,d such that for all v : ^ R, the action of f over v is given by 


f{v):= {fo,v)N+ 




d 


= /o(a:/iV)w(a;/iV) +^ fk{x/N)d^^v{x/N) W{{x + Ck/n)-W{x) 


^ / J \ ' ]^d 

a:6T^ 

Moreover, = ||/o|||r + Y.k=i (T^)' denote the functional f by 


k=l x£T}j 

II X 112 


f = k-jZ€jk- 

k=l 

Proof, let u be the unique weak solution obtained by Lemma 14.41 when A = 1 , A the identidy matrix 
and, / £ In particular, we have (/, u) = {u,v)i^w,N■ Thus, it suffices to consider fo = u and 

fk = d^^u, where fc = 1 ,..., d. □ 

4.2. Connections bet-ween the discrete and continuous Sobolev spaces. Let us consider the natural 
partition of induced by T^. For each we denote by Qn{x) the set 

Qn(.x) = {j/ £ T'^; Xk <yk < Xk + 1/N, k = 1,2,. .. ,d} . 

Let un ■ —>■ R be a mesh function. We will now introduce some interpolation different schemes, that 

is, procedures to extend un from to the continuous torus T'*. The first scheme is the piecewise-constant 
interpolation given by 

UN{y) = X! UN{x)l{QN{x)}{y). 

^^T7'^N 

However, extend a function un ■ R to in such a way that it belongs to Hi is not 

straightforward. To do so, we need the definition of W-interpolation, which we give below. 

Using the standard construction of d-dimensional linear interpolation, see for instance [9], it is possible 
to define the W-interpolation of a function un '. jf^N Rj with W(x) = Wk{xk) as defined in ID- 

More precisely, let y £ Qn{,x), then the W-interpolation of it^v, namely, is given by 


ulkiy) = un{x)-\-'^ d^,^UN{x){Wk{yk)-W{xk)) - 

fc=i 

+ ■ • ■ ^w,un{.x) W{W,{y,) - W^{xi)) 

k —1 i —1 

i^k 

d 

+ • -d^^UN^x) Y[0^k{yk) - Wk{Xk)). 

fc=l 


Let US consider a new type of interpolation, which are actually d (one for each m = 1,..., d) interpolations: 





d 

u^N^y) = un{x) + '^ d^^UN{x){Wk{yk)-W{xk)) + ■ ■ ■ 

k^l 

k^m 

d 

+ dwi ' " 9{^rn-i^Wm+i ' " ^Wd'^Nix) JJ ((y^) — (x^)). 

i^l 

i^m 

A simple calculation shows that 

(15) ^ (y)- 

We now establish the connection between the discrete and continuous spaces by showing how a sequence 
of functions defined in can converge to a function defined in T'^. 

We will say that the sequence of mesh functions (/at )Nen in or (T^), k = 1,... ,d, converges 

strongly (weakly) to the function / € L^(T'^) or L^^(T'^), as TV —>■ oo, if /at —>■ / in L^(T‘^) or L^^(T'^) 
strongly (weakly), respectively. Similarly, a sequence (uAr)ArGN in Hi converges strongly (weakly) to 

u £ if —>■ M strongly (weakly) in Hi,wi^‘^)- 

With respect to convergence of functionals, we say that a sequence of functionals /at £ with 

representation 

d 

In = fo,N - ^ d^^fk,N 

fe=i 

converge weakly (resp. strongly) to / £ if the sequence of functionals given by 

d 

Fn = fo,N - ^ dxkfk,N 
fe=l 

converges weakly (resp. strongly) to / in i/^^(T'^). 

Lemma 4.8. Let Uff : be a sequence of functions such that: 

• There exists a constant C > 0 satisfying 

(16) l|wAf|li2(T^) = ^ N>1, 

xeT^ 

and, one of the sequences {u’^), (un) or (u^^) (for some m = I,... ,d) converges weakly in L'^{T^) 
to a function u : —>■ R os TV —> oo. Then, the others also converge to u in the same manner. 

• For j £ {1,2,..., d}, there exists a constant Cj > 0 satisfying 

(17) I“a^IIl^.(t^) = ^Nix/N)[W{{x + ej)/N) - W{x/N)] < Cj for od > 1 

^ xGT^ 

and, one of the sequences (itj^), (un) or (u^^) (for some m = 1,... ,d) converges weakly in (T'^) 
to a function u : R. as N ^ oo. Then, the others also converge to u in the same manner. 

Proof. We start with the L^(T‘^) case. Suppose (uj^) converges weakly to u in L^(T'^). We will prove that 
{un) also converges to u. The remaining cases can be handled in a similar manner. Note that (1161) implies 
the uniform boundedness of the norms ||Mj^||A2(T'i)j ll^wII l2(T‘^) IIl 2(T‘^) because the values of the 

functions u'^, un and in each cell Qn{x), x £ T^, lies between the largest and smallest values of UAf 
at the vertices of Qn{x). Hence, each of theses sequences is weakly compact in L^(T^). From this, to prove 
the convergence we are interested, it is enough to prove the convergence un^ dy —>■ fj,j dy for all 
4) £ (^““(T^) (notice that we have, by assumption, that dy —>• fj.a u4> dy). 
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Let <I>Ar be the piecewise-constant function which coincides with $ at the lattice-points ;^T^. By uniform 
convergence of ‘Pat to we have dy —■ fj^ u<!> dy. Let us consider 

Rn ■■= / {uli-UN)^Ndy= ^ d>(x) / {u% -un) dy. 

1 T<i JQn(x) 

By explicit form of the interpolations, we have 
Rn = 




^ $(a:){ ^a^^MAr(x) / {Wk{yk) - W{xk)) dy + ■ ■ ■ 
fc=i 

d p d 

'^d^,---dw;,^^d^^+,---dw^uN{x) / Y[{Wi{y^) - W^{xi)) dy 

1 _1 Qwix) •_i 


'Qn{x) 

i^k 


P d 

+ d^r---dw^uN{x) Y[{W^{y^) - Wi{x^)) dy] . 

d Qr^(x) .•_i ^ 


Qn{x) 2—1 

After a summation by parts we transfer in each term inside the braces one of the LL-differences from un 
to $ and thus, the resulting expression is 

Rn = — 


where 


N = - [ (Wkiyk)-Wixk))dy + --- 

i-=l JQn(.x) 

<L f^i + W I 

^\_ /(#)-/(^) 


i=l ■ 


gN 


/(v) 


-, for X G T^, 


and the exchange of the order of integration is due to Fubini’s Theorem. Note that 

^ l - Xi .+ l/N 2 

/ [Wi^ ivij ) - Wi- {Xi^ )]dy^j < 

j = l dx,. 


^Nix) 

i=i' 

■ "^Wi,^N{x) YllWi^ixy^ + 1/N) - Wi.(x2j]) 

j=i 

1 2 ^ 

$(x)[LF 22 (X22 + 1/iV) - Wy, (Xy, )]) Un(x) njVF., (X2, + 1/iV) - W., (Xy, )])' 


3=2 


i=l 


where in the last equality the constant C > 0 does not depend on TV. 

In this way, \Rn\ is bounded by a finite sum (that depends on d) of terms of the form 


NN^ 




1 / 2 , 




2=1 


^ '^N 
1/2/ 1 


2=1 


1/2 








where the previous estimates follow from Cauchy-Schwarz inequality. By hypothesis (|16l) and the fact that 
$ G C'°°(T'^) the sums inside the brackets are bounded. Thus, Rn goes to zero as TV —>■ oo. Therefore, we 
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have shown that 


un^n dx 


w<i> dx. 


/fd 


jfd 


Hence 


/ MAr$ dx= ' 

JTl Jjii 




Un^ dx = / UN^Ndx+ / UN[^ — ^N]dx^ / dx. 


/jd 


This concludes the proof for the first case. Now, let us consider the (T'^) case. The proof is similar to 
previous case. To keep notation simple, fix, without loss of generality, j = 1. The beginning of the proof is 
exactly the same as in the previous case, one just needs to replace the Lebesgue measure dy by the product 
measure d{yi <S> Wi). The essential difference comes in the estimation of the term 

r ^ 2 

(d^ ^ix)d^.^ ■ ■ ■ un{x) / [W,^ {y,^) - (xi ^)] d{yi ® bhi)) . 

^ '^QNix) ^ 

We will now provide this estimation. The previous term is bounded above by 


1 ^ 

+ 1/iV) - TTi(xi)] + 1/N) - lT.,(x.,)])' 

i=i 

= (+ ^) _ X 

X {dwi, ■ ■ ■ W7V (x) P [Wi^ {xi. + ^) - Wi, (x.,)]) 

3=2 

^ f[Wl{xi + j^)-Wi{xi)] ^ 1^2^ ^ ^2 

<(- {Ci}_^u^{x + e,/N)j . 


So the new term |i?Ar| is, again, bounded by a finite sum (that depends on d) of terms of the form 

1 ^ , 1, _n, ,Nl/2/ 1 __ 1, , ,,,,,^1/2 

N 






N 


By hypothesis (HZD and the fact that $ € C'°°(T‘^) the sums inside the brackets are bounded and thus Rn 
goes to zero as iV ^ 00 . Therefore, we have shown that 


Td 


UN^Nd{yi(SiWi)^ / (i(2/i (g) Wi). 


Ijd 


Hence 


/fd 


UN^ d{yi ^Wi) = / un^n d{yi (^Wi) 


jjd 


I'd 


un[^ - $Ar] d{yi ® Wi) 


jfd 


d{yi ®Wi). 


This concludes the proof. 


□ 


We will now state and prove the strong counterpart to the above Lemma: 

Lemma 4.9. Let : T'^ —>■ K. fee a sequence of functions and C > 0 be a constant such that 

(18) ll«Jvlli/i,w(T^) for all N >1. 

Then, if one of the sequences (un) or[u^^'^) (for some m = 1,... ,d) converges mL^(T‘^) to a function 

u : as N ^ 00 , then the others also converge to u in the same manner. 

Proof. Let us consider the case u% —>■ m in L^(T'^) and we will prove that —> u in L^(T‘^) . The remaining 

cases can be proved analogously. Recall the notation introduced at the beginning of this subsection and note 
that 


(19) 


Rn = f {uj^ - dy < y] f {u^! - dy . 

1 w JQn{x) 
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A simple computation show that the quantity \u*pf{y) — u^^\y)\ assumes its largest value on Q{x) at one 
of the vertices of Qn{x) and that this value is equal to {Wm (Xm + l/N) - Wm{xm))\d^l^UN\- That is, 

( 20 ) 'maXy^Q^(^^)\u*j^{y) - u^^\y)\ = + l/N) - Wm(x^))'maXy.\y_^\=iiM, y„=x,n\^Wm'^N{y)V 

In fact, 

d 

Uwiy) - U^N\y) =(Wm{ym) - Wm{Xm)^ ^d(;(^^UN{x) + ^ d(;(,.d^^UN{x){Wj{yj)-W{Xj)) H- 

j=X,j^m 

d d 

+ ^ ^Wi ■ ■ ■ ^Wj-i^Wj+i ■ ■ ■ 9w^d^^uN{x) {Wi{yi) — Wi{xi)) 

i=l 


a 

''' ^Wd^Wm.'^Nix) {Wi{yi) — Wi{xi))'^ = 

2=1; i^m 

( 2 /m ) (^m)^ (Z/l j * ' * ; ym — 11 2 /m+l; ‘ ; Ud) 


where (d^^UN)* denotes the W— interpolation of the function O^^un in the {d — l)-dimensional setup. 
In this way, i?v is bound above by 




[ ([Wm{xm + ^/N)-Wjn{xm)]'max^.\^_^\^i/N^^^^x^\d^^UN{z)\\ dy = 
-^(\Wm{Xm + l/N)-WmiXm)]maX^.\^_^\^l/N,z^=x^\dw^UNiz)\^ < 

N 

c Wm{l)-Wm{0) (^dN^UNiz)y[WmiZm + l/N)-Wrn{z)]j^ = 




Cl 


W^{l)-W,n{0),^N „2 ^ Wm{l) - Wm{0) .. ,,2 


N 


'\\^Wm'^N\\wm,N — Cl 


N 


\\un\\h. 


1,W(T^) 


where in the previous expression, Ci is a constant that depends on d (one may take, for instance, Ci = 2'^). 
So, by (dH), limAT^oo Rn = 0. 

Thus, it follows that 

f (it^^ -ufdy <2 f {u% - + {ulf - uf dy 0, 

dfi JTi 


as iV —>■ oo. This concludes the proof. 


□ 


Finally, we have the following Lemma regarding strong and weak compactness: 

Lemma 4.10. Let un : T'^ —>■ M 6e a sequence of functions such that there exists some constant C > 0 such 
that, for every N > 1, 

(21) jv(T^) < C. 

Then, the sequence (u^) forms a uniformly bounded set in iJi ^/(T’^), which is therefore strongly precompact 
in L‘^(T‘^) and weakly precompact in Hi^w0^‘^) ■ 

Proof. To prove the lemma, it is enough to show that (1211) implies the uniform boundedness of M/(T<i) 

and then we apply m Proposition 2.9]. 
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Note that 


|wAfllL 2 (T-i) = / iu*j^iy)fdy= f {u%{y)f dy 


< 


E 


\z-x\=i/N 


Y {un{z))^ - 7^ E = 2‘^l|uAr||i2(T^) 


zeTij 


where z, in the previous expression, represents the summation over all vertices of the cell Qn{x). Further, 
for m = 1, 2,..., d, we have, by an explicit computation. 


9w„u%{y) =dY'>J'N{x)+ Y ^Wj^Wm'^N{x){Wj{yj) -W{xj)) -\ - 

d d 

+ E + 1 n 

j=l-,j^m i=l 

d 

+ dw,--- dw^dYuN{x) (Wi(?/i) - W,{xi)) = 

2=1; iz^m 

(yij*** ; I/m— 1 5 ym + 1; * * ‘ ; 

where denotes the hF— interpolation of the function O^^un in the (d — l)-dimensional setup. 

In this way, 

\\dw^ul,\\Y iT^) = [ (dw^uUy))^ d{y^ ^Wrn) = 

Y f < Y [ E idY^N{z))^d{y^®W^) 

2 X ^ / ^AT / \ \ 9 {/ Xm “(“In 


< 


7Vd-i 

Thus we have shown that 






where the last inequality follows from (EH), and this concludes the proof. 


□ 


We will now obtain the converse procedure, that is, how to use a measurable function / in L^(T'^) to 
properly define a mesh function in T^. This is not straightforward, since the restriction of / to the set 
{0, 1/N,... ,{N — 1)/7V}'^ is not well-defined (this is a set of Lebesgue measure zero). 

For arbitrary functions / £ L^(T‘^) and g € consider the mesh functions defined on the 

discrete torus of the form 

( 22 ) fN{x) =N^ [ f{y)dy, 

d Qn(x) 

and 

+ UN)'- ® 

where x € and Qjv(cc) = {y G T^; 0 < i/i — Xi < N~^, i = 1,... ,d}. 
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For a generalized function f € in the form (O, / = /g — J2k=i ^xkfk, where /g G L^(T‘^), and 

fk € we consider 

d 

(24) fNix) = fo,N{x) - d^Jk,N{x), 

k^l 

where x G ;^T^, /g^^r is the discretization given by (|2^ . and /fe_Ar {k = 1,... ,d) is the discretization given 
by (1^ . /at clearly belongs to 

The following lemma shows that this procedure yields a good approximation for /. 

Lemma 4.11. Let f G and g G (T*^) for some j = 1,2,... ,d. Let Jn and gM be the mesh 

functions defined in equations (j22p and (j23l) . Let fN and'gN be the piecewise-constant interpolation for f^ 
and gpf. Then, 

11/ ~ /Af||L2(T‘i) 0, 

and 

lls ~ 9 n \\ l ^ . (T'i) 0, 

as N ^ oo. Ln particular, if F G then F/v —>■ F strongly in 


Proof. We will prove the first assertion. The proof of second is analogous. Recall the notation used in the 
previous lemmas. Since the continuous functions are dense in ( and also in LF^w-O^'^) ) it is enough 

to consider / continuous. Thus, 


dy = 


II/-/iv|lL2(Td) = / (fiy)-fNiy)) dy = Y i (fiy)-N‘^ [ f{z)dz\ 

Y[ (iV"* / [f{y)-f{z)]dz) dy < Y i f if iv) - f dz dy , 

f.^^jQNix'^) ^ JQnIx'^) ' tC 1 T‘i JQn(x)JQn{x) 

N "iV 

where the last inequality follows from Holder’s inequality. Thus, the previous expression is bounded above 

by 

dv f [f{y + v)-f{y)?dy<2‘^ Y / if (v + v) - f iv)? dy . 

JQn(x) JQn(x) 


Y 




,<i/Ar 




To conclude, note that by compacity of T'’* the continuous function / is, in fact, uniformly continuous. So, 
for a fixed e > 0, there exists Nq such that ||?7||oo < l/-^o implies that \f{y + g) — f{y)\ < e^/^/2'’*. Therefore, 
the previous expression is bounded by e, and the proof of the Lemma follows. □ 


5. Homogenization 

Our main goal in this Section is to prove the convergence of energies, namely Proposition 15.61 and also 
the Homogenization Theorem, Theorem 15.81 This last result is presented with fairly general hypotheses on 
the matrices . In Proposition 15.91 we provide an example of a very large class of functions that admit 
homogenization and in subsection 15.31 we consider a scenario in which represents a random environment. 
We begin with some definitions and auxiliary results. 

5.1. Definitions and auxiliary resnlts. We now focus on the analysis of the asymptotic behavior of the 
sequence (un) given by solutions of the equations 

(25) \UN - V^A^V^un = f^, A > 0 

where f^ are fixed functions defined on ;^T^, = {d^^,... ,8^ = {d^^,... ,d^^) are the 

difference operators, and A^ = {akk)dxd are diagonal matrices, A^ = {akk)dxd, of order d satisfying the 
ellipticity condition: there exists a constant 6 > 0 such that 6~^ < a^(a;) < 9, for every a; G and 
k = 1,... ,d. 
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The continuous counterpart of the theory developed in Subsection |TT] can be found in m- More precisely, 
one can find results on existence, uniqueness and boundedness of weak solution of the problem 

(26) Auq — = /, A > 0. 

We say that the diagonal matrix = (a^) i/-converges to the diagonal matrix A = (akk), denoted by 
A^ A, if for every sequence of functionals on Hi and / € such that —>■ / as 

iV —)• oo strongly in we have 

• un —t uq weakly in Hi^w (T'^) as fV —>■ oo, 

• ^kk^Wk'^N -t akkdwkUQ weakly in (T'^) for each /c = 1,..., d, 

where um : —>• R is the solution of the (1251) and uq € Hi^iy(T'^) is the unique weak solution of the (1261) . 

In this case, we say that the diagonal matrix A is a homogenization of the sequence of random matrices A^. 
We also say that the operator \7AS/\y is a homogenization of the sequence of random operators A^V^. 
Next, we present an estimate that will be useful in this Section. 


Lemma 5.1. Let f as N ^ oo weakly in Then, if un : ^ R is the solution of the 

(1^ . there exists C > 0, not depending on N, such that 


jv(T^) < C. 

Proof. Using lemma l4^ we obtain a constant Ci, not depending on N such that 


Note that fN f weakly in means that fN ^ f weakly in Since ||/iv|| \ = 

||/||^-i(.j,d), and /at is weakly convergent, it is bounded, and thus there exists a constant C 2 > 0, such that 


Thus, for all iV > 1, 


ll^^llffi,w(T^) — ^lC '2 = C. 


□ 


We will now prove a very simple version of the compensated compactness Theorem. 

Lemma 5.2 (Compensated Compactness in (T'^)). For k = 1,... ,d, let {gk,N) md (wk^N) be two 

seguences of funetions on (T'^) such that 

9k.N ^ gk, strongly in and Wk,N ^ Wk, weakly in 

where gk,Wk e Then, 

gk.NWk,N gkWk weakly in 


Proof. Let (f S S)vv(T'^), then 


Ifd 


gk,NWk,N(fd{x^ ® Wk) = / {gk,N - gk)wk,N 4 >d{x^ ® Wk) + / gkWk,N<pd{x'^ ® Wk). 


Ifd 


Ifd 


Let us deal with each term in the right-hand side of the previous equation. 
Note that 


Ifd 


{gk.N - gk))wk,N(l)d{x^ ® Wk) 


< \\gk,N - gk\ 




Note that, from Lemma 13.51 cj) is bounded, and since {wk,N) is a weakly convergent sequence, its norm is 
uniformly bounded. Therefore, the previous equation tends to zero as fV —>■ 00 . 

We now deal with the other term. Begin by recalling that </> is bounded. Thus, qk4> € since 

is a set of finite d{x'^ (8> IUfc)-measure. Therefore, the weak convergence of Wk,N to Wk implies that 


Ijd 


gkWk,N(l)d{x^ ®Wk) = / Wk,N{gkf')d{x’^ ®Wk) 
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WkPkf'dix^ ® Wk). 




This concludes the proof. 


□ 


We conclude this subsection with a version of Compensated Compactness Theorem for discrete approxi¬ 
mations. 

Corollary 5.3 (Compensated Compactness for Discrete Approximations). For k = 1,... ,d, let {qk,N) and 
{vk,N) be sequences of functions on such that 

qk,N ^ qk, strongly in and Vk,N ^ Vk, weakly in ilk 

where qk,Vk € Then, 

qk,NVk,N -t qkVk weakly in (T'^). 

Proof Let Wk,N = Vk,N and gk,N = qk,N- Thus, from definition, Wk,N converges weakly in (T'*) to 

Vk, and gk,N converges strongly in to qk- Therefore, we may apply lemmato conclude that 

Vkj^,N = Vk,Nqk.N = Wk,Ngk,N Vkqk weakly in 

□ 

Remark 5.4. One should notice that lemma 1 5. ill is, indeed, a version of the Compensated Compactness 
Theorem. In fact, the classical assumptions would be gk,N gk and Wk,N Wk weakly in (T'^), 

and dx,^gk,N —^ h strongly in where dx,^gk,N should be understood as in lemma 1X71 However, 

in our setup, the functional induced by gk,N coincides with —dxkPk.N, and thus, since \\dxgk,N\\H~^(jd'j = 
\\gk atIIls CTd-i, we recover strong convergence for gk n- The reason why the functional induced by gk n is 
~dxkgk,N instead of the standard functional is because the standard functional given by 

gk,N<l>dx 

Jf<>- 

is not well-defined for gk,N G fact, gk,N may not belong to L^(T'^). 

5.2. Main Results. We are now in a position to state and prove the homogenization of the difference 
operators introduced in the previous Section. We begin by proving an auxiliary lemma that will be needed 
in the proof of the convergence of energies of a sequence of homogenized matrices. 

Lemma 5.5. Let X be a Banach space and X* its dual. If fN € X* is such that fiq —>■ / € X*, and un G X 
is such that un ^ u G X weakly. Then, 

fN{uN) f{u), 

as N ^ oo. 


Proof. Since un converges weakly in X, it forms a bounded sequence, that is, there exists C > 0 such that, 
for each N >1, HuAfU < C. We also have 

\fN{uN) - f{uN)\ < II/at - f\\x*\\uN\\x < C\\fN - f\\x*, 

which tends to zero as iV —>■ oo. On the other hand, from weak convergence, /(uat) —>■ f{u) as iV —> oo. 
Therefore, fxiux) —t f{u) as —>■ oo. □ 


The next proposition shows that even though the iJ-convergence only requires weak convergence in its 
definition, it yields a convergence in a strong sense (convergence in the L^-norm for the piecewise-constant 
interpolation). 


Proposition 5.6. Let A, as N ^ oo, with un being the solution of (I25p . where f € H^/(T'^) is 

fixed, f^ ^ f strongly in and, uq is the weak solution of (EH). Then, the following limit relations 

hold true: 

Mat —t Mo in L^(T'^), 


1 


ceTi, 


u%{x/N) 


N—¥oo 


Jjd 


UQ{x)dx, and 
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— X! X! ^kk{^l^){dwk'^N{x/N)Y [Wk{{xk + l)/iV) - Wk{xk/N)] / afefe(x)(avi^^Mo(a;))^rf(a^''®W4). 


Nd 

Proof. We begin by proving that 
(27) 


fc = l ' 


Jd 


f^{uN) /(■Uo), 


as —> cx). As we plan on using Lemma |4.81 we need to obtain a bound. By Lemma l5.11 the sequence un 
is II • ||i,vv bounded uniformly. In particular, there exist constants Ci, > 0, /c = 1,..., d, such that, for all 


iV > 1, we have 


(28) 

and 

(29) 


^ y: «i(i/jv) < Cl, 




Nd 


— idw,UN{x/N))^{Wk{{xk + 1)/N) - Wk{xt/N)) < 


ceTt 


Now, observe that from Lemma HTTl there exist functions /o,Ar,..., fd,N such that 

d 

fN{uN) = (/o,Af,MAf)L2(x^) + y~!(/fc L'^ (Td)- 


fe=l 


This motivates us to define the following functionals gi^N G ^w^C^w) by go,N{v) = (/o.at,'k)z, 2 (t^), and 
9k,N{v) = {fk,N, d^ v)i ^2 „d ), for fc = 1,..., d. Note that 

“ x>=<S)W^' 

go,N{v) = {fo,N,v)mTXy and gk,N{v) = (fk,N ,d^v) ^2 (T-i)- 

We have, by hypothesis, that /at —>■ / strongly, which means 

d 

II/a( — /ll^^i(Td) = ll/Af,0 — /o|lL2(T<i) + ^ \\fk,N — (T-^) b, 


k=l 


as N ^oo. Thus, go,N -t /o strongly in ^^(T'^), gk,N fk strongly in (T"*). 

Now, observe that un —t uq weakly in H^{Td) means that u’^ —t uq weakly in id^(T‘^), which in turn 
implies that 

(30) u% —>• Uo, weakly in L^(T'^), 

and 

dwkU*N ^ dwkUo, weakly in 
Nevertheless, by equation (US, we have 

dWkU*N = {dwk'd-N)^^^ , 

and thus 

(31) ^ dwkUo, weakly in (T"*). 

Using Lemma IT^ (1^ implies 

(32) UN —t Uo, weakly in 

and (1^ implies 

(33) d^^UN^dwuUo, weakly in 

Now, since 50 ,tv -t /o strongly in and gu.N -t fk strongly in we may apply Lemma 


5.51 to equations (13^ and (1331) to obtain that 
(34) 30 .Ar(uAr) —>■ /o(uo). 
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as iV —>■ 00 and. 




(35) 

Note that 


9k,N{d^^UN) fk{dwkUo), as iV -)• oo. 

d 

f{uo) = /o(mo) + '^fk{dwk{uo)), and 


k=l 


fN{uN) = go,N{uN) + y^,9k,N{dw^{uN)) = go,Af(uAr) + gk,N (d^^UN)- 

k^l k^l 

Thus, from (f34|l and (f3^ . fN{uN) converges to /(no). Now, note that 

/(mo) = A / uldx + y^ akkidwkUo)^d{x'"®Wk), 


H 


since, by the hypothesis that A —i> A, ug is the weak solution of Auq — VAVwUo = /■ Note also that 
/^(uat) = (XuNjx/N) -V^A^V^un{x/N))un{x/N) 


ceT^ 


A 

w 


<{x/N) - T ^ un{x/N)V^A^V^un{x/NI 

ceT^ x(^T% 

which, after a summation by parts in the above expressions, and using that fN{uN) /(uq), we obtain 


^ ^ u%{x/N) + ^ ^ ^ a^kidw,^N{x/N)nWki{xk + 1)/N) - W^Xk/N)] 


Nd 


ceTl 


(36) 


k=lxGT% 
N—¥C !0 


A / uldx + 'V] / akkidwkUofdix'^'SiWk). 

JTi JT<i 


Suppose that un does not converge to uq in L^(T'^). That is, there exist e > 0 and a subsequence (unj) 
such that 

IImaTj — Mo||L 2 (Td) > e, 

for all j. By Lemma [4. 101 we have that there exists v G L^(T‘^) and a further subsequence (also denoted by 
ujv„) such that 

5(5 3 ^9P • r2/nr'd\ 

—> X, inL (T ). 

Using Lemma 14.91 we further obtain that 

UNj —^ n, in L (T ). 

This implies that 

UNj —>■ weakly in L^(T‘^), 

but this is a contradiction. Indeed, from H convergence of Ajv to A, we have that un ^ uq weakly in 
iJi vvCIT'^), which means that u’^ —>■ Uq in iJi^vv(T‘^)- Thus u*j^ —)• Uq weakly in L^(T‘^). Using Lemma 14.81 
this implies that uat —> uq weakly in L^(T'^). In particular, 

UNj —i” Uq, weakly in L^(T'^), 

and at the same time ||z; — uo||/, 2 C]['<i) > £■ Therefore, un —> uq in L^(T'^). The proof thus follows from 
expression (1551) . 

Corollary 5.7. Let uq G (Du/(T'^) and un '■ j^'^n ® such that limuAr = uq in L^(T'^). Then, 

||mo — >■ 0 . 


□ 
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Proof. We have that un —t uq in L^(T'^) means that 

(37) IIUAT - Mo||L2(T<i) -5-0, 

a.s N ^ oo. On the other hand, from the definition of the set Qm{x) in Subsection 14.21 and, the fact that 
the functions in are right-continuous (see Remark I3.4L we have that 


Mo It^ 


Uo 


pointwise as fV —>■ oo. From Lemma 13.51 wq is bounded, and thus uq is bounded, and therefore integrable 
(T"^ has finite Lebesgue measure). Thus, we can use the Dominated Convergence Theorem to conclude that 


(38) ||mo It^ — Mo||L2(T<i) -1 0, 

as iV —>• oo. Therefore, using (|37)) and (|38l) . we obtain 

II^Ar — = llwiv — Mo It^ l|L2(T‘i) 

< IImat — Mo||l2(T<1) + I|mo — Mo It^ ||L2(Td) -5-0, 

as N goes to oo. □ 


We will now state and prove the main result of this paper. 


Theorem 5.8. Let = {a^^)k=i:d be a sequence of diagonal matrices and 9 > 0, such that 6 ^ < 9, 

and 1/a^ converges weakly in Then, admits a homogenization. 

Proof of Theorem \5.8[ Fix / G and consider the problem 

(39) XuN - V^A^V^un = fN, 

where /at is the discretization obtained in equation From Lemma 14.111 fN ^ f strongly in Pdf^lff^). 

Using Lemma l5.11 there exists a unique weak solution un of the problem above such that its i7)\y-norm 
is uniformly bounded in TV. That is, there exists a constant C > 0 such that 

From Lemma 14.101 there exists a convergent subsequence of un (which we will also denote by un) such that 

UN —t M, weakly in iJi^vv(T'^). 


In particular, 


(40) d^^UN^-^dw^u weakly in Ll,^^^{T‘^). 

Applying (15^ to un, we obtain 

d 

A||ma/||a2(t^) + ^ = fN{uN) < \\fN\\H^^(T'f)\WN\\Hl^,{^%)■ 

k=l "" 

Thus, for each k = 1,... ,d, and using Lemma 15.11 we have 

< WfNWn^ATf) ■C'- 

Using the same argument (from weak convergence of the functionals) we can find a constant Ci > 0 such 
that < Cl. Therefore, 

hkk^Wk'^NWL^ (T^) < CiC. 

This, in turn, implies that 

W^kkdw.UNWLl^^^^iT^) = \\akkdw,UN\\Ll^^^^(Tf^) 
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is uniformly bounded in N. Thus, since is a separable Hilbert space, we can find a further 

subsequence (also denoted by un), such that 

(41) ^ vo,k weakly in 

as iV —>• oo, uo,fc being some function in (T'^). 

Since (ujv) is uniformly bounded in the Sobolev-norm and L^(T‘^) is precompact in this space we have 

UN ^ u strongly in T^(T‘*) . 

In particular, 

un ^ u strongly in 

On the other hand, {\un — A^V^un) converges strongly (to /) in Therefore, 

strongly in if 

From the very definition of the functional A^V^un, the previous convergence means that for each k, 

akkdw^UN vo^k strongly in 

Denote by bk € the weak limit of the sequence (1/a^): 

l/ttkk bk weakly in 

Since 9~^ < a^(a;) < 9, we have that is uniformly bounded, and > 9. Further, bk > 0 

d{x^ 0 kFfe)-a.e. . In fact, if there exists a measurable set A with positive d{x^ 0 Wk) measure, such that 
= 0 in A, take the function (/)= 1a- Then, 

0 < 0measure(H) < J —(()(i(x^ 0 H4) —>■ J 6fclAd(a;^ 0 II4) = 0. 

A contradiction. Thus bk > 0. 

From the Compensated Compactness Theorem (take = a^kb^Wu'^N and = I/a^ in Corollary 15.31) : 

-^Ukkdwk'^N bkVo,k, weakly in 
“fcfc 

On the other hand, 

-^a^k^Wk'^N = d^^UN -)■ dw^u, weakly in ■ 

“fefc 

From uniqueness of the weak limit, we have that = bkVo^k- Since bk ^ 0, we have that 

vo-k = ^dw^u- 
Ok 

Thus, we can summarize our findings: 

UN Strongly in (T'^) , 

dwkUN ^ dwkU weakly in and 

Ok dwkUN strongly in (T'^). 

Therefore, u solves the problem 

Xu — VAVwm = /, 

where A is the diagonal matrix with entries given by 1/bk- 

To conclude the proof it remains to be shown that we can pass from the subsequence to the sequence. 
This follows from uniqueness of weak solutions of the problem (1251) , see m Proposition 3.4]. The fact that 
any converging subsequence is a solution to the same problem, and the fact that ujv is uniformly bounded 
in the Sobolev norm, thus we can find a convergent subsequence (thus a sequence that do not converge to 
the solution, must converge to somewhere else, since the limit point is also a solution, uniqueness shows the 
result). □ 
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We will now provide an example of a very large class of functions that admit homogenization. Recall the 
definition of the space Mw given below Remark l3.II 

Proposition 5.9. Let A = (akk) G Mw be a diagonal matrix such that 0~^ < Ukk for some 6 > 0 and, the 
discretization = (a^) be the sequence of diagonal matrix obtained from (1231) . Then, the sequence 
admits a homogenization. 

Proof. It is clear that 0 < 0“^ < a^. Further, it is clear that the right-continuity implies the pointwise 
convergence of to Okk- Finally, from Lemma [3751 the functions Okk are bounded, and thus the sequences 
are bounded. From the Dominated Convergence Theorem 

akk akk strongly in 

On the other hand, we also have the pointwise convergence of 1/a^ to Ifokk, and the bound 1/a^ < d 
implies that 

1/a^fc 1/afcfc strongly (T"')- 

Therefore, the result follows from Theorem l5.8l □ 

5.3. Homogenization of Random difference operators. In this subsection we consider the homoge¬ 
nization problem when the matrix A^ represents a random environment. More precisely, we focus on the 
analysis of the asymptotic behavior of the sequence (un) given by solutions of the equations 

XuN - V^A^V^un = f^, 

where f^ are fixed functions defined on ;^T^, ... ,d^^) and ,...,5^^) are the 

difference operators and, the random diagonal matrix A^ = {akk)dxd of order d represent the statistically 
homogeneous rapidly oscillating coefficients. Therefore these equations are driven by the random difference 
operators A^V^, and to fully understand them, we need to understand the random matrices A^. Thus, 
let /i) be a standard probability space and {T^ : —>• 0;a; € be a group of J^-measurable and 

ergodic transformations which preserve the measure 

• Tx : Lt ^ LI is J^-measurable for all x 

• = ^(A), for any A G J" and x G 

• To = J, Ta; o Ty = Ta^+y, 

• Any / G L^{Ll) such that /(Tata) = /(w) /r-a.e., for each x G Z^, is equal to a constant /i-a.e.. 

Note that the last condition implies that the group is ergodic. We call the underlying probability space 
random environment, and a point a; G D a realization of the random environment. 

Let us now introduce the vector-valued J^-measurable functions {6fefe(a;);fc = l,...,d} such that there 
exists 0 > 0 with 

(42) 0-1 < bkkitd) < 9, 

for all a; G D and k = 1,... ,d. Then, define the random diagonal matrices whose elements are given by 

(43) b^^{x) -.= bkkiT^xO}) , X G , k = l,...,d. 

Let us show some weak convergences associated the random environment (0^) defined in (H5|) . First, note 
that by Birkhoff Ergodic Theorem, we have, 

(44) b’^f. — E[bkk\ weakly in L^(T"*) a.s., 

for fc = I,..., d. We will need an similar result for (T*^). 

Denote by measure induced by function Wk. By Lebesgue decomposition, there exist, function 

g such that, 

TWk = gX + 

where gX and A"*" are singular measures and A denotes de Lebesgue measure. Let 14 C T be the support of 
A-*- and = Tx...xTxVfexTx...xTc 14 in the A:-th component. 
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Define af : 


as 


(45) 


-‘fcfc 


(X) = 


6^fc(x) if 2J'= n giv(x) = 0, 

E[bkk] if 03'= n giv(*) 0 

Here Qn{x) denotes the cell in T”* with vertex in a; G 

Qn{x) = {y & <yi- Xi < l/N} . 

Define the random diagonal matrices = (a^) and consider the problem 
(46) \UN - 

The main result of this subsection is 

Theorem 5.10. Let A^ he a sequence random matriees, as defined previously in dUl). Then, almost surely, 
A^iuj) admits a homogenization, where the homogenized matrix A does not depend on the realization to. 

Proof. The proof follows from Lemma [5.111 that ensures the hypothesis of the Theorem l5.8l are valid for this 
sequence. □ 

We conclude this subsection with following result. 

be as defined above. Then, 

E[hkk\ weakly in a.s., 


Lemma 5.11. Let 

N 

^kk ~ 


and 


l/affc —^ Bkk weakly in a.s., 

where the function Bkk is given by 

B U) = */03'=ngAr(x) =0, 

\l/E[hkk\ */Q3'=ngAr(x) ^0. 

Proof. Following the notation introduced above, let 23^ = UQj'iinQN(a:)#0 Qn{x). We have, 




pointwise. 


Therefore, 

(47) 


1®^ I®'' 


strongly in T^(T='), 


_^ 

where 53 stands for the closure of the set 53'=. Let </> G SwCIT'^) be fixed. By Lebesgue decomposition we 
have, 

[ a^k^'dix^ ®Wk) = [ a^i.(j)gdx + f a^(()(i(a;'= ® A'*'). 

JT<i Jt‘‘- 

Note that the support of the measure d(a;'= (8) A-*-) is conhned in the set 53'=, defined above. Since 53'= C 53^, 

we have that is almost everywhere constant, namely = E[bkk], with respect to the measure d{x^^X-^). 
Thus, the second integral in the right-hand side in previous expression is equal to 

[ E[bkk](l>d{x'^ ® X-^) . 

JTd 

For other side, by and (H71) . we have 
(48) —y Elbkkjl^^ky weakly in a.s., 

here X'^ denotes the complementary set of X. 

So, the first integral in the right-hand side in above expression is equal to 


/jd 


(^kk^<s%(l^9dx 


y^^kk^[<a%Y^9dx = J^^E[bkk]l<aY^gdx + = 




From previous convergence, (ITTl) and (H5)) . the right hand-side in previous expression converges to 

E[bkk]^m’‘^9dx+ E[bkk]^r^kY<l}gdx = / E[bkk](t>gdoi 

JTd. L J 


tX . 


Jfd 

Then, we have showed that 


lim / 

N —yoo J'jpd 


® Wk) = 


Ifd 


E[bkk]4'd{x'' <S) A-*-) 


Ifd 


E[bkk]4>9dx = / E[bkk]4'd{x'' <S) Wk) 


iTd 


and this concludes the proof of the first statement. 

The second statement follows directly from the first by observing that 1 /bkk is also an ergodic sequence, 
and everything we did above may also be done to the sequence 1 / . □ 

6. Application 

To conclude the paper we will provide an application of a new result on probability theory which is an 
improvement of the result obtained in [14] in two directions: first, it considers a more general model; second, 
it has a more natural and simpler proof. It is also noteworthy that the homogenization results obtained in this 
paper were the key results in proving the main Theorem in the article [2] . For each choice of diagonal matrix 
function satisfying the hypotheses of Theorem 15.81 we have a corresponding version of hydrodynamic 
limit. For instance, we may obtain a version in random environment using Theorem l5.101 and also a version 
with any diagonal matrix function A € 


6.1. The hydrodynamic limit. We will now use the standard vocabulary in probability theory, that 
is, cadlag functions means right-continuous functions with left limits; tightness is a property regarding 
compactness; weak convergence is actually weak* convergence. The reader is also referred to [ 6 ] and references 
therein. 

We will begin by recalling some definitions. Recall in ([TJ the definition of function W, consider a sequence 
of operators that satisfies the hypothesis of the Theorem l5.8l We will consider, for instance, the 

random environment introduced in subsection 15.31 For each x € and j = 1,..., d, define the symmetric 
rate ^x,x+ej ~ ^x+ej^x by 


af{x/N) af{x/N) 


N[W{{x + e,)/N)-W{x/N)] N[W,iix,+1)/N)-W,ix,/N)y 
where af{x) is given by (1451) . and ei,..., Cd is the canonical basis of Also, let b> — 1/2 and 

Cx,x+ej (r?) = 1 -I- b{T]{x - ej) -|- r]{x + 2ej)} , 

where all sums are modulo N. 

Distribute particles on in such a way that each site of is occupied at most by one particle. Denote 
by T] the configuration of the state space { 0 , 1 }^" = {77 : —>• { 0 , 1 }} so that r]{x) = 0 if site x is vacant, 

and r]{x) = 1 if site x is occupied. 

The exclusion process with conductances is a continuous-time Markov process {rjt : t > 0} with state 
space {0, 1 }^", whose generator acts on functions / : {0, 1 }’’'''^ —> K as 


(50) 


{LNf){r]) = X! X] 


d:,x+ej \ _ 


V) - fiv)} ) 


where is the configuration obtained from rj by exchanging the variables ri(x) and r](x + ej): 

ri{x + Cj) ii y = X, 


(51) 




V{x) 

My) 


a y = X + Cj, 
otherwise. 


We consider the Markov process {rjt : t > 0} on the configurations {0,1}^^ associated to the generator 
Lpf in the diffusive scale, i.e., Ljv is speeded up by N'^. 

We now describe the stochastic evolution of the process. Let x = {xi,... ,Xd) & T^. After a time given 
by an exponential distribution, at rate ^x,x+ejCx,x+ej{r]) the occupation variables r]{x) and r]{x + ej) are 
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exchanged. Note that only nearest neighbor jumps are allowed. If W is differentiable at x/N € [0,1)'^, 
the rate at which particles are exchanged is of order 1 for each direction, but if some Wj is discontinuous 
at Xj/N , it no longer holds. In fact, assume, to fix ideas, that Wj is discontinuous at Xj/N, and smooth 
on the segments {xj/N,Xj/N + eej) and {xj/N — eej,XjfN). Assume, also, that Wk is differentiable in a 
neighborhood of Xk/N for k ^ j. In this case, the rate at which particles jump over the bonds {y — ej,y}-, 
with yj = Xj, is of order 1/A^, whereas in a neighborhood of size N of these bonds, particles jump at rate 
1. Thus, note that a particle at site y — Cj jumps to y at rate 1/N and jumps at rate 1 to each one of the 
2d—l other options. Particles, therefore, tend to avoid the bonds {y — ej,y}. However, since time will be 
scaled diffusively, and since on a time interval of length N'^ a particle spends a time of order N at each site 
y, particles will be able to cross the slower bond {y — ej,y}. The scaling limits of this interacting particle 
systems in inhomogeneous media may, for instance, model diffusions in which permeable membranes, at the 
points of discontinuities of the conductances W, tend to reflect particles, creating space discontinuities in 
the density profiles. For more details see [14]. 

The effect of the factor Cx,x+ej iv) is the following: if the parameter b is positive, the presence of particles 
in the neighboring sites of the bond {a:, x + ej} speeds up the exchange rate by a factor of order one, and if 
the parameter b is negative, the presence of particles in the neighboring sites slows down the exchange rate 
also by a factor of order one. More details are given in Remark 16.21 below. 

Let A = {ajj)dxd be a diagonal matrix belonging to Mu/ with ajj > 0,j = l,...,(i, and recall from 
subsection [31 the operator defined on J)u/(T'^): 

d 

VAVw '■= dxjajjdwj, 

1=1 

where A S Mu/. 

A sequence of probability measures {yiN : > 1} on {0,1}’'^^ is said to be associated to a profile 

po : [0, 1] if for every 5 > 0 and every function H € S)vv(T'^): 


(52) 


lim upf 

N—¥(yD 




^ y] H{x/N)ij{x)- f H{u)po{u)di 


>6 


= 0 . 


Let 7 : T'^ ^ [Z, r] be a bounded density profile and consider the parabolic differential equation 


r dtp = VAVwHp) 
\ P(0,-) = 7(-) 


where the function 4> : [/, r] —>• ffi. has bounded derivative, its derivative is also away from zero, and t G [0, T], 
for T > 0 fixed. 

A function p : [0, T] x —)• [1, r] is said to be a weak solution of the parabolic differential equation (1551) if 
the following conditions hold. $(p(-, •)) and p(-, •) belong to L^([0, T], iLi_u"(T'^)), and we have the integral 
identity 


/id 


p{t,u)H{u)du — / p{0,u)H{u)du = 




I Id 


^{p{s,u))VAVwH{u)duds 


for every function H G 2)u^(T‘^) and all t G [0,T]. 

Existence of such weak solutions follow from the tightness of the process proved in subsection 16.31 and 
from the energy estimate given in [141 Lemma 6.2]. Uniquenesses of weak solutions was proved in [12] . 

The main result of this Section is the following. 


Theorem 6.1. Fix a continuous initial profile po : —>■ [0,1] and consider a sequence of probability 

measures pN on {0, associated to po, in the sense of (l52l) . Then, for any t > 0, 


lim P 

N —¥oo 





H(u)p{t, u) du 



= 0 


for every S > 0 and every function H G (Du^(T‘^). Here, p is the unique weak solution of the non-linear 
equation (|53p with / = 0, r = I, 7 = po and <I>(q;) = a + aa^. 
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Remark 6.2. The specific form of the rates Cx,x+ei is not important, but two conditions must be fulfilled. 
The rates must be strictly positive, they may not depend on the occupation variables ri{x), r]{x + ei), but they 
have to be chosen in such a way that the resulting process is gradient (cf. Chapter 7 m [B] for the definition 
of gradient processes). 

We may define rates Cx^x+e, to obtain any polynomial of the form $(«) = a + 'Yli 2 <j<m ^ m > 1, 
with 1 + '^ 2 <j<md ^3 ^ Z®*' instance, m = 3. Then the rates 

Cx.x+efip) = Cx,x+ei{'n) + b {ri{x - 2ei)ri{x - cf) + pix - ei)r]{x + 2ei) + ^{x + 2ei)r]{x + ici)} , 

satisfy the above three conditions, where Cx,x+ei is the rate defined at the beginning of Section 2 and a, b are 
such that 1 + 2a + 36 > 0. An elementary computation shows that <i>(a) = 1 + aa^ + ba^. 

6.2. Proof of Theorem 16.11 A simple computation shows that the Bernoulli product measures {i/^ : 0 < 
a < 1} are invariant, in fact reversible, for the dynamics. The measure is obtained by placing a particle 
at each site, independently from the other sites, with probability a. Thus, is a product measure over 
{0,1}'"''^ with marginals given by 

: p(x) = 1} = a 

for X in T^. For more details see [6l chapter 2]. 


Consider the random walk {Xt}t>o of a particle in induced by the generator given as follows. Let 
f,x,x+ej be given by (l49l) . If the particle is on a site x G T^, it will jump to x + Cj with rate N‘^fx,x+ej- 
Furthermore, only nearest neighbor jumps are allowed. The generator L^r of the random walk {X(}t>o acts 
on functions / : —>• R as 


where. 






/(- 


d 

E 

j=i 


‘-W(v) 


)-/( 


H" ^x—e 




It is not difficult to see that the following equality holds: 


(54) 


L^/(x/iV) = 


1=1 


iafd^Jfix) 


:= S/^A^V^fix). 


The counting measure m^v on iV is reversible for this process. This random walk plays an important 
role in the proof of the hydrodynamic limit of the process pt- 

Let II(R+, {0,1}^«) be the path space of cadlag trajectories with values in {0,1}’"''^. For a measure p,N 
on {0,1}^'^, denote by the probability measure on II(M+, {0,1}^'^) induced by the initial state and 
the Markov process {pt ■ t > 0}. Expectation with respect to P^^ is denoted by 

Let A4 be the space of positive measures on with total mass bounded by one endowed with the weak 
topology. Recall that nf G Ai stands for the empirical measure at time t. This is the measure on 
obtained by rescaling space by N and by assigning mass 1 /N‘^ to each particle: 

(55) = 7^ E St{x)Sx/N , 

where (5„ is the Dirac measure concentrated on u. 

For a function H : T‘^ ^ M., {ir^, H) stands for the integral of H with respect to : 

{'xf.H) = ^ H{x/N)pt{x) . 


This notation is not to be mistaken with the inner product in L^(T‘^) introduced earlier. Also, when ttj has 
a density p, Tr{t,du) = p{t,u)du. 

Fix T > 0 and let D{[0,T],M) be the space of A4-valued cadlag trajectories tt : [0,T] —>■ Af endowed 
with the uniform topology. Note that A4 is endowed with the weak topology, which is metrizable, since 
is a compact metric space. Thus, this uniform topology is well-defined. For each probability measure pN on 
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{0,1}^", denote by the measure on the path space £)([0,T], A^) induced by the measure and the 

process introduced in (1551) . 

Fix a continuous profile po : T‘^ ^ [0,1] and consider a sequence {pN : iV > 1} of measures on {0, 
associated to po in the sense (l52l) . Further, we denote by the probability measure on £>([0,T],Ad) 
concentrated on the deterministic path 7r(t, du) = p{t, u)du, where p is the unique weak solution of (1531) with 
7 = Po: Ik =0, Tfc = 1, k = 1,..., d and $(a) = a + ba^. 

In subsection l6.3l we show that the sequence : TV > 1} is tight, and in subsection l6.4l we characterize 

the limit points of this sequence. 


6.3. Tightness. The goal of this subsection is to prove tightness of sequence : TV > 1}. Fix A > 0 

and consider, initially, the auxiliary Al-valued Markov process {Ilf’^ : T > 0} defined by 

xeZ'^ 


for H in 1)^(1"^), where iVf is the unique weak solution in iJi ,w{%) of 

Aidf - = \H- VAVwH, 

with the right-hand side being understood as the restriction of the function to the lattice T^, which is 
well-defined, since H G Dyy{T^) and from Remarkwe have that VAVwH belongs to Thus 

the right-hand side belongs to C'w(T‘^) and it is right-continuous. Remark 13.41 

We first prove tightness of the process : 0 < f < T}. Then we show that ■ 0 <t <T} and 

{tt^ ■ 0 < t < T} are not far apart. 

It is well known [6] that to prove tightness of {Ilf’^ : 0 < f < T} it is enough to show tightness of the 
real-valued processes {Ilf(id) : 0 < t < T} for a set of test functions ii : —)• R dense in C'(T‘^) for the 
uniform topology, for instance we can use 2)vv(T‘*). 

Fix a function ii : T'^ —R in (Dw(T'^). Keep in mind that nf’^(ii) = and denote by 

the martingale defined by 

(56) TWf’^ = - [ dsN^LNiTT^.H^) . 

Jo 

Clearly, tightness of nf’'^(ii) follows from tightness of the martingale TV/f^’^ and tightness of the additive 
functional J* ds N'^LN{Tr^, H^). 

A long computation, albeit simple, shows that the quadratic variation of the martingale M^A 

is given by: 

E E [d^,H^{x/N)f[W{{x + e,)/N) - Wix/N)] x 

i=l aieT'' 

In particular, by Lemma |4.4l 


/ Cx,x+ej{Ps) [ns{x + ej) - Ps{x)]'^ ds . 
Jo 




< 


Nd 


El 

i=i 


H 


N\\2 


< 


\ \\Wj,N — 


C{H)t 
XNd ’ 


for some finite constant C{H), which depends only on H. Thus, by Doob inequality, for every A > 0, 5 > 0, 


(57) 


lim P, 
N —¥oo 




sup 

0<i<T 


M, 


N,\\ 


> s 


= 0 


In particular, the sequence of martingales {M^A ; TV > 1} is tight for the uniform topology. 

It remains to examine the additive functional of the decomposition (1561) . The generator of the exclusion 
process Ljv can be decomposed in terms of the generator of the random walk Ljv. By a long but simple 
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computation, we obtain that N'^Liy{TT^, H^) is equal to 

^^NHx)ix/N)Tj{x) - (^NHx)ix/N){Tooh2,j)iri)} 

J = 1 xGT% xGT^ 

xeT% 


where {tx '■ x G Z'^} is the group of translations, so that {Tx'rj){y) = r]{x + y) for x, y in Z*^, and the sum is 
understood modulo N. Also, /iij , /i 2 j are the cylinder functions 


KAv) = vAAAj) , h2,j{y) = y{-ej)y{ej) . 

Since is the weak solution of the discrete equation, we have by Remark 14.51 that it is also a strong 
solution. Then, we may replace by = A^x ~ + VAVwH in the previous formula. In 

particular, for all 0 < s < t < T, 


/ 


drN^LN{K,HA) 


< 


(l + 3|6|)(t-s) 

Nd 




\UA{x/N)\ 


It follows from the estimate given in Lemma 14.41 Remark 13.21 Lemma 13.51 and from Schwartz inequality, 
that the right hand side of the previous expression is bounded above by C{H,b){t — s) uniformly in N, 
where C{H,b) is a finite constant depending only on b and H. This proves that the additive part of 
the decomposition (1551) is tight for the uniform topology and, therefore, that the sequence of processes 
{nf’^ : > 1} is tight. 


Lemma 6.3. 


The sequence of measures 


: > 1} is tight for the uniform topology. 


Proof. Fix A > 0. It is enough to show that for every function H G (Dw(T'^) and every e > 0, we have 


lim ¥ N sup I - {ttA, H)\> e 

v-s-oo Lo<t<r 


= 0 , 


whence tightness of follows from tightness of 11^’^. By Chebyshev’s inequality, the last expression is 
bounded above by 



UN 


sup 

.0<t<T 



H\\ 


2 

L2(T^)’ 


since there exists at most one particle per site. By Theorem 15.81 Proposition 15.61 and Corollary 15.71 WH^ — 
iL||^ 2 (.j,rf ^ 0 as A^ —oo, and the proof follows. □ 


The reader should compare the proof given in this Section to the proof given in |14) , to check that the one 
given here follows closely the standard approach given, for instance, in [5], whereas the approach given in 
|I4) is non-standard. This means that the theory provided in this article made the study of hydrodynamic 
behavior of exclusion processes with conductances tractable in the standard fashion. Thus, future work on 
this field will be simplified. 


6.4. Uniqueness of limit points. We prove in this subsection that all limit points Q* of the sequence 
concentrated on absolutely continuous trajectories 7r(t, du) = pit, u)du, whose density p(t, u) is a 
weak solution of the hydrodynamic equation (I53p with / = 0, r = I and 4 >(q:) = a -|- aa^. 

We now state a result necessary to prove the uniqueness of limit points. Let, for a local function g : 
{0,1}* ^ M, g : [0,1] —>■ R be the expected value of g under the stationary states: 

gia) = [gig)] . 
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For £ > 1 and d-dimensional integer x = (xi,... ,Xd), denote by ri^{x) the empirical density of particles in 
the boxB^(x) = {{yi,...,yd) £ ;0 < yt - Xi < £}: 

= ji Y. ■ 

yGB^(x) 

Let Q* be a limit point of the sequence and assume, without loss of generality, that converges 

to Q*. 

Since there is at most one particle per site, it is clear that Q* is concentrated on trajectories iTtidu) which 
are absolutely continuous with respect to the Lebesgue measure, 7rt(du) = p{t,u)du, and whose density p is 
non-negative and bounded by 1. The reader is referred to li Chapter 4] for further details. 

Fix a function H £ 2)vv(T‘^) and A > 0. Recall the definition of the martingale Mj. introduced in the 
previous Section. From (1571) we have, for every i5 > 0, 


lim P„ 
N—¥oo 


sup > S 

0<t<T 


= 0 , 


and from (1551) . for fixed 0 < t < T and i5 > 0, we have 


lim 


N —¥oo 


- «,iFf) - fdsN^L^{7r^,H^) 
Jo 


> S 


= 0. 


Note that the expression N‘^Ln{tt^,H^) has been computed in the previous subsection in terms of 
generator L^r. On the other hand, — XH VAVwH- Since there is at most one particle 

per site, we may apply Theorem 15.81 along with Proposition 15.61 and Corollary 15.71 to replace ,H^) and 
{ttq,H^) by {nt,H) and {tto,H), respectively, and replace by VAS/wH plus a term that vanishes 

as N ^ oo. 

Since E,y^[hij] = o?, i = 1, 2 and j = 1,... ,(i, we have by the replacement Lemma [TH Corollary 5.4] 
that, for every t > 0, A > 0, i5 > 0, i = 1, 2, 

rt 


lim lim sup Pu^ 

e^O 


Y '^AVwH{x/N) X 


> <5 


= 0 . 


Since (x) = e {Yl‘i^i[xj/N,Xj/N -|-eej]), we obtain, from the previous considerations, that 


lim lim sup Q 


£—^0 


N—><X) 


W,N 

fJ-N 


pt ^ 

“'o i=i 


> 5 


= 0 . 


Using the fact that converges in the uniform topology to Q*, we have that 


lim ( 
£—^0 


{nt,GxH) - {t:o,GxH)- [ ds (^{e ‘^7r5(TT[-, •-f ecj])) , Ua) 

do \ ! 


> (5 


= 0 . 


Recall that Q* is concentrated on absolutely continuous paths TTt(du) = p{t,u)du with positive density 
bounded by 1. Therefore, £~’^tTs (ntih- -I- ecj]) converges in L^(T^) to p(s,.) as e j, 0. Thus, 


> (5 


= 0 . 


{TTt.H) - {tto,H) - [ dsmps),VAVwH) 

Jo 

Letting i5 j, 0, we see that, Q* a.s., 

f p{t,u)H{u)du — f p{0,u)H{u)du = f f ^{p{s,u))'\7AX/wH{u)duds . 

JT'^ Jt^ Jo Jt^ 

This identity can be extended to a countable set of times t. Taking this set to be dense we obtain, by 
continuity of the trajectories ttj, that it holds for all 0 < t < T. 

From |141 Lemma 6.2], which we may easily adapt to our setup by using the uniform ellipticity condition 
(1421) of the environment, we may conclude that all limit points have, almost surely, finite energy, and therefore. 
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by [m Lemma 4.1], •)) G ?"], iJi_iv(T‘^)). Analogously, it is possible to show that p(-, •) has finite 

energy and hence it belongs to L^([0,T],iLi_^^(T^)). 

Proposition 6.4. As N oo, the sequence of probability measures converges in the uniform topology 

to Qiv. 

Proof. In the previous subsection, we showed that the sequence of probability measures is tight for 

the uniform topology. Moreover, we just proved that all limit points of this sequence are concentrated on 
weak solutions of the parabolic equation (1531) . The proposition now follows from the uniqueness proved in 

m- □ 

Proof of Theorem \6.1[ Since converges in the uniform topology to Qw, a- measure which is concen¬ 
trated on a deterministic path, for each 0 < t < T and each continuous function iL : R, 

converges in probability to Jj^dup(t,u)H{u), where p is the unique weak solution of (l53l) with Ik = 0, 
rk = 1, 7 = po and d>(Q:) = a + aa^. □ 


7. APPENDIX 


In this appendix we will present, for the reader’s convenience, part of the proof of the [S] Lemma 1]. This 
result was used in Proposition 13.31 


Proposition 7.1. The space 2)w(T) is dense in C'(T), the space of continuous functions in T, in the sup 
norm, || • ||oo- 

Proof. Let / : T —R be a continuous function, and e > 0. From uniform continuity, there exists 5 > 0 such 
that \f{y) — f{x)\ < e whenever \x — y\ < S. Choose an integer n > 5~^ and consider the function p : T —>■ M 
defined by 

/([j + l]/n) -/(j'/Af) 


9{x) = 


Wk{[j + l]/n) - Wkij/N) 


if + l)/n]}{x) , 


with 1{A} being the indicator function of the set A. g can be seen as a discrete derivative of / with 
respect to Wk- Thus, it is natural to guess that integrating this function with respect to Wk would yield 
an approximation of /. Thus, let G : T —>• R be given by G{x) = /(O) -I- ,,,^g{y)Wk{dy). By the very 

definition of g, G{j/n) = f{j/n) for 0 < j < n. 

Since n > S~^, we can use the definition of G to obtain, for jfn < x < (j -I- l)/n, 

\G{x)-f{x)\ < |G(a;) - G(j7n)| -h |/(a;) -/07n)| < 2e , 
whence ||G — /||oo < 2e. Note that 


(58) 


'( 0 . 1 ] 


gdWk = 0 . 


It remains to be shown that the function G may be approximated in the sup norm by functions in 
Note that the only restriction we had when choosing the set {0,1/n,..., (n — I)/n} is that the distance 
between two consecutive points is less than S. Therefore, we may replace these n points, by any such points 
satisfying this restriction. Since Wk is strictly increasing, it has a countable number of discontinuities, and 
then, we may assume, without loss of generality, that Wk is continuous at the points {0,1/n,..., (n — I)/n} 
(in the sense, that we may replace these points if it is not the case). Let {Hm : to > 1} be a sequence of 
smooth functions Hm : T —>• R, with \Hm{x)\ < UpHoo, for all x € T, and such that \iu].mHm{x) = g{x) for 
xn 7 2). Then, the Dominated Convergence Theorem implies 


(59) 


lim f \Hmiy) - g{y)\dWkiy) = 0. 

m—^oo Jr^ 


Let {Fm : TO > 1} be the sequence of functions Fm '■ T 


defined by 


Fm{x) = /(O) / H{y)dWk{y). 


'(O.a:] 
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Therefore, it is immediate that 


||G-F„|U< [ \H^{y)-g{y)\dWkiy)'"-^0. 

Jt 


This shows that we can approximate, in the sup norm, any function / G G(T) by functions of the form of 
Fjn- To conclude that is dense in C'(T), we must show that To this end, note that 

Fm{x) = /(O) + + J H'^(z)dz'^dWk{y) 

= /(o) + bmWk{x) + [ dWk{y) [ H'^{z)dz , 

J(0.x] Jo 


where bm = ^fm(O) — Wfe(l) ^ Hm{y) dWk{y)- Since Fl^ is continuous, FI'^ G Gvi/(T). Then, using the 
characterization of one may conclude that Fm belongs to TiWkj for each m > 1. □ 
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